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Abstract. The relativistic hydrodynamic model is applied to describe the expansion of the dense 
matter formed in relativistic heavy-ion collisions. The hydrodynamic expansion of the fluid, supple- 
mented with the statistical emission of hadrons at freeze-out gives a satisfactory description of the 
observables for particles emitted with soft momenta. Experimental data for the transverse momen- 
tum spectra, elliptic flow and interferometry radii give constraints on the properties of the fluid, its 
equation of state and viscosity coefficients. The role of the fluctuations of the initial profile of the 
energy density is discussed. 
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INTRODUCTION 

Experimental studies of heavy-ion collisions at ultrarelativistic energies have the aim at 
creating a drop of very dense and hot matter and measuring its properties. Measurements 
at the BNL Relativistic Heavy Ion Collider (RHIC) cover a range of energies up to 
^/sNN = 200GeV for Cu-Cu and Au-Au systems. Recent results from the CERN Large 
Hadron Collider (LHC) extend the available energies to ^/spj^ = 2.76TeV for Pb-Pb 
interactions. 

The results accumulated in the last decade indicate that a fireball of strongly inter- 
acting, expanding fluid is formed [1, 2, 3, 4, 5] in the collision. Observables related to 
the emission of particles with soft momenta demonstrate the onset of a strong collective 
flow of matter. Particle spectra are written using Fourier coefficients of the directed vi, 
elliptic V2, triangular V3, or higher harmonic flow 

+V3COs(3(,(l-'i'3))+...) ■ (1) 

The flow coefficients can be extracted using the orientation angles ^ or from many- 
particle cummulants. Transverse momentum spectra of identified particles can be quan- 
titatively understood as coming from the convolution of the collective flow of the fluid 
elements and of the thermal emission [6, 7]. This observation is even more striking for 
the azimuthal asymmetry of momentum distributions, the elliptic flow. The observation 
of the elliptic flow has a very convincing interpretation in terms of the collective ex- 
pansion of a source with an azimuthally asymmetric initial geometry [8]. Two-particle 
correlations of same particles allow to measure the interferometry radii of the emission 



region [9] . The values of the radii and their dependence on the momentum of the pion 
pair show the presence of a strong transverse flow. The extracted short emission time 
points towards a rapid expansion of the fireball. 

We discuss the interpretation of the expansion of the fireball in terms of the perfect 
fluid [10, 11, 12, 13, 14, 15, 16] and viscous hydrodynamics [17, 18, 19, 20, 21, 22]. 
Many efforts have been devoted to the construction of a model giving a realistic descrip- 
tion of the data. These studies yield a dynamical picture of the space-time evolution 
of the bulk of the matter created in the collision. At the same time they give access to 
some physical properties of the system, the size and life-time of the fireball, the equa- 
tion of state, the shear viscosity coefficient. Many unknowns that come into play in the 
construction of the model, require the use of as many experimental constraints as possi- 
ble to limit the possible systematic errors. In that way, a satisfactory description of the 
soft observables limits the effect of the model uncertainties or even helps to get some 
additional information, e.g on the role of the fluctuations in the initial state. 



RELATIVISTIC HYDRODYNAMICS 

The expansion of the fireball is described using relativistic hydrodynamics. The basic ob- 
ject in the formalism of second-order viscous hydrodynamics is the energy-momentum 
tensor [23] 

^ ^ py^^v _ pg^v ^ + nA'^^ (2) 

the energy-momentum tensor of the perfect fluid gets a correction in the form of the 
stress tensor Tt^^ from shear viscosity and from bulk viscosity n, is the fluid velocity, 
= giiv _ j^Mj^v noiQ, that additional transport coefficients are possible in general 
[23, 24]. The hydrodynamic equations 

9^7^" = (3) 

and the equation of state connecting the pressure p and the energy density e are supple- 
mented with dynamical equations for the stress corrections 

A-AV..H.., = '-n^r_^ - (^) (4) 



urd,n = Z^z^^i^ _ ^n^5, ( ^ ) . (5) 



Tn 2 Tn V CT^ 

Two viscosity coefficients appear in the equations, the shear viscosity 7] and the bulk 
viscosity ^ and two relaxation times Tji and Tn- The value of the viscosity coefficients 
can be estimated in kinetic theory if the cross sections are known. However, for very 
strongly interacting systems kinetic theory breaks down [25]. In AdS/CFT calculations 
the ratio of the shear viscosity to entropy reaches a very small value 77/5 = l/4;r 
[26]. The small value of the shear viscosity, as inferred from the experimental data 
[1, 2, 3, 4, 5], justifies the use of the perfect fluid hydrodynamics as a first approximation. 
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FIGURE 1. Transverse momentum spectra of pions produced in Au-Au collisions at = 200GeV at 
different rapidities, compared to BRAHMS Collaboration Data [33] (from [16]) 



Perfect fluid calculations can explain a large number of observations [10, 11, 12, 13, 
14, 15, 16, 27]. Calculations have been performed both in simplified, boost-invariant 
2+ 1-D geometry as well as in 3 + 1-D. The generation of the transverse flow is cor- 
rectly described as the acceleration of the collective flow from pressure gradients in the 
source. The predicted elliptic flow overshoots the experimental points, the calculation 
can be made closer to the data if a hadronic cascade stage is introduced after the hy- 
drodynamic expansion [27]. The Hanbury Brown-Twiss (HBT) correlations could not 
be described in early calculations using an equation of state with a first order phase 
transition [12, 28, 29, 30]. The agreement is greatly improved when using lattice QCD 
inspired parameterization of the equation of state with a cross over [14, 31]. 

To illustrate the quality of the predictions of perfect fluid hydrodynamics , we show an 
example of a 3 + 1-D hydrodynamic simulation [16] with a cross-over equation of state 
[32]. In Fig. 1 are shown results for the pion spectra at different rapidities. Spectra of 
pions and kaons can be well described for soft momenta p± < 1 .5GeV in a perfect fluid 
3 -|- 1-D model. Up to now, there exist only one calculation of viscous hydrodynamics in 
3 + 1-D [22], other viscous codes use the boost-invariant geometry. 

We illustrate the sensitivity of the HBT radii, and especially of the ratio Rout /Rside to 
the softening of the equation of state. The increase of this ratio has been proposed as a 
measure of the strength of the first order transition [34]. Experiments at RHIC and the 
LHC exclude such a behavior. When using the equation of state with minimal softening 
around Tc (Fig. 2), one obtains the best agreement with the measured interferometry radii 
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FIGURE 2. Temperature dependence of the square of the velocity of sound, for a cross-over equation 
of state (sohd hne) [32] and for an equation of state with a softening around Tc (dashed hne) (from [16]) 
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FIGURE 3. HBT radii for hydrodynamic calculation with a cross-over transition (dashed lines) and for 
an equation of state with some softening around Tc (solid lines) (from [16]). 



(Fig. 3). The success of the simulations using a cross-over equation of state [14, 31] is 
a phenomenological confirmation of recent lattice QCD calculations of the equation of 
state [35]. 

The directed flow has been measured at ultrarelativistic energies. The formation of the 
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FIGURE 4. Directed flow of charged particles from 3 + 1 -D hydrodynamic calculation, compared to the 
data of STAR and PHOBOS CoUaborations [36, 37]. (from [38])- The color band indicate the systematic 
error form the range of the uncertainty in the initial tilt of the source. 



flow in the hydrodynamic expansion implies a non-boost invariant geometry and requires 
a mechanism breaking the symmetry around the beam axis. Asymmetric emission from 
forward and backward going participants [39] leads to a tilt of the source. The model 
reproduces the observed value of the directed flow, in Au-Au and Cu-Cu collisions (Fig. 
4). The formation of the directed flow can serve as probe of the early non-equilibrium 
imbalance between the longitudinal and transverse pressures [40]. 

The elliptic flow is generated in the expansion of an initial azimuthally deformed 
source. The final elliptic flow depends on the time of the evolution, on the viscosity 
coefficient [41] and on the initial eccentricity. The measured elliptic flow can be used 
to extract the value of the shear viscosity coefficient [42, 43]. The main uncertainty is 
related to the value of the initial eccentricity of the fireball. Depending on the assumed 
model of the initial density profile (Glauber model or KLN model) the extracted value 
of the shear viscosity lies between rj/s = 0.08 and rj/s = 0.2. The calculation of the 
elliptic flow for identified particles requires a realistic description of the hadronic stage 
of the collision [44, 45]. If this stage is modelled by hydrodynamics, bulk viscosity must 
be introduced. In Fig. 5 the results for V2 of different particles are shown using ri/s = 0.1 
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FIGURE 5. Elliptic flow of identified particles scaled by the initial eccentricity in hydrodynamic 
calculations compared to the data of PHENIX and STAR Collaborations [46], [47], [48] (from [44]). 




FIGURE 6. Elliptic flow coefficient in Au-Au coflisions at = 200GeV (STAR data [48]), and for 
Pb-Pb collisions at 2.76TeV (ALICE data [5]), compared to viscous hydrodynamics results [49]. 



and C,ls = 0.04 in the hadronic phase. In summary, hydrodynamic studies of RHIC data 
have shown that a strongly interacting matter is formed in the collision, with a cross-over 
equation of state and a small viscosity. The expansion of the the fireball is very rapid. 



FROM RHIC TO LHC 

From the recent data on Pb-Pb collisions at the LHC [50, 5] a similar picture of the 
dynamics emerges as at RHIC energies. The multiplicity and the life-time of the system 
increases, but the global properties are very similar. In a first approach, the same param- 
eters as used at ^J~s = 200GeV may be used (rj/s = 0.08 for Glauber initial conditions 
and C,/s = 0.04 below Tc) increasing only the initial energy density. Several calculation 
exist for LHC energies [51, 52, 53, 45, 49, 44, 54, 55]. The elliptic flow of charged par- 
ticles can be well reproduced (Fig. 6). The spectra of identified particles compared to 
calculations are shown in Fig. 7. The agreement is satisfactory. However, for light parti- 
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FIGURE 7. Transverse momentum spectra of identified particles from viscous hydrodynamics, com- 
pared to preliminary ALICE Collaboration data [56]. 
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FIGURE 8. HBT radii for pions emitted in Pb-Pb Collisions at = 2.76GeV. The hydrodynamic 
model results are denoted by the dashed lines for zero initial flow and by the solid lines when the initial 
pre-equilibrium flow is included [44], ALICE Collaboration data [50] (from [44]) 



FIGURE 9. Distribution of wounded nucleons in a Glauber Monte-Carlo event, the orientation of the 
participant plane is rotated with respect to the reaction plane, (from [61]) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

[GeV/c] 



FIGURE 10. Elliptic flow of ;r+, and p in Pb-Pb collisions, from hydrodynamic calculations 
compared to ALICE Collaboration data [62] . 



cles the mean transverse momenta seem too large. The HBT radii can be accounted for 
by the viscous hydrodynamic calculation (Fig. 8) with the transverse flow generated in 
the expansion of the fluid. It shows, that the size and the life-time of the system is well 
understood in dynamical models. 

The shape of the initial distribution in the transverse plane (Fig. 9) is largely influ- 
enced by fluctuations [57, 58, 59]. As a measure of viscosity effects the triangular flow 
is especially important [60] . In particular, it is difficult for the model to fit at the same 
time the elliptic and the triangular flow. Another problem, is to reproduce the mass split- 
ting in the dependence of V2 between pions and protons. The elliptic flow of identified 
particles is reasonably well described (Fig. 10), except for more central events [62]. In 
Fig. 11 we illustrate the difficulties of the model to reproduce both the strength and the 
mass splitting of the triangular flow. Some hints of the problems with common predic- 
tions for the spectra of pions, kaons and protons is already visible in Fig. 7. It indicates 
that the very strong transverse collective flow seen at LHC is not entirely accounted 
for in present hydrodynamic calculations. As noted before, the spectra and the ellip- 
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FIGURE 11. Same as Fig. 10 but for triangular flow. 



tic flow of identified particles are sensitive also to the final hadronic rescattering. More 
advanced approaches to fluctuating initial conditions use event by event hydrodynamic 
calculations [57, 22]. 

In summary, the hydrodynamic expansion is a realistic description of the behavior 
of the bulk matter created in relativistic heavy-ion collisions. Viscous hydrodynamics 
reproduces the observed transverse momentum spectra, the elliptic flow of charged 
particles and the HBT radii. One can conclude that the hot matter created in the collision 
behaves like an almost perfect fluid with a cross-over like equation of state. Fluctuations 
of the shape of the initial state are evidenced in the azimuthal asymmetry of the emitted 
particles. Some discrepancies remain in the description of the momentum dependence 
of the elliptic an triangular flows of identified particles. 
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